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Abstract 

Digital signal processing techniques were employed to investigate the joint use of charge division and risetime analyses for 
the resistive anode (RA) coupled to a microchannel plate detector (MCP). In contrast to the typical approach of using 
the relative charge at each corner of the RA, this joint approach results in a significantly improved position resolution. 
A conventional charge division analysis utilizing analog signal processing provides a position measured resolution of 170 
^m (FWHM). By using the correlation between risetime and position we were able to obtain a measured resolution of 
92 p,m (FWHM), corresponding to an intrinsic resolution of 64 ^m (FMHM) for a single Z-stack MCP detector. 
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1. Introduction 

Position-sensitive microchannel plate (MCP) detectors 
are a powerful and widely used tool in imaging of electrons, 
photons and ions [l|. Since their inception in the late 
1950’s 0, they have been used in a variety of applications 
ranging from ion-molecule reaction scattering experiments 
0 to fast neutron radiography Q. In this detector, an 
incident electron, photon, or ion ejects an electron from 
the micorchannel plate. The microchannel plate acts as 
a continuous dynode electron multiplier with millions of 
independent channels providing a typical amplification of 
sslO^. Two MCP plates (chevron configuration) or three 
plates (Z-stack configuration) can be stacked to achieve 
higher gains. Different methods exist for measuring the 
position of the electron cloud exiting the microchannel 
plate thus determining the position of the incident par¬ 
ticle. The princ^al techniques are : multi-anode [^, he¬ 
lical delay line induced signal 0, and resistive anode 
0 . The resistive anode (RA) technique is particularly ap¬ 
pealing for its simplicity. The electron cloud emanating 
from the MCP stack is incident on a two-dimensional re¬ 
sistive sheet. Readout of the charge at the four corners 
of the sheet provides a measure of position of the inci¬ 
dent particle via charge centroiding. Utilizing this simple 
approach a position resolution of 134 p,m (FWHM) for 
a chevron conhguration has been achieved [9|. Although 
for many application a chevron stack of two microchan¬ 
nel plates provides sufficient amplification, the detection 
of low-intensity signals near the one electron limit requires 
configurations involving three (Z-stack) or more MCPs. 
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For a Z-stack configuration the best reported resolution is 
100-200 fim (FWHM) 0- By use of more complex MCP 
stack configurations and through the use of retarding po¬ 
tentials, position resolution as good as 50 ^m (FWHM) 
has been realized [niEp. In this work, we investigate 
the use of pulse shape analysis to improve the position res¬ 
olution obtained with a simple Z-stack MCP-RA detector. 


2. Experimental setup 


Depicted in Fig. [T] is the experimental setup used 
to determine the position resolution of the MCP-RA de¬ 
tector. Alpha particles from a ^'^^Am radioactive source 
(A) are emitted towards a secondary electron emission foil 
(C). Electrons are emitted by the passage of an a particle 
through the 1.5 /rm thick aluminized mylar foil. Ejected 
electrons are accelerated to an energy of 1.3 keV by the 
potential difference between the aluminized foil and a wire 
grid (D). The electrons then pass through slits in a stain¬ 
less steel mask (E) directly in front of the microchannel 
plate stack (F) before impinging on the front surface of 
the MCP-RA detector. The precision mask, fabricated by 
laser micro-machining 1^, has ten slits each measuring 
100 /im by 7620 p,m and spaced by either 4.2 or 4.5 mm 
apart. Electrons that pass through slits in the mask are 
amplified by the Z-stack MCP and detected on a resistive 
anode (G). Alpha particles traverse the aluminized mylar 
foil and are detected by a fast scintillator/photomultiplier 
tube (PMT) assembly (B) placed directly behind the alu¬ 
minized foil. The MCP used was a standard Z-stack (APD 
3 40/12/10/12 D 60:1) with 10 iim diameter microchan¬ 
nels provided by Photonis USA [l^ which was coupled to 
a 40 mm diameter RA from Quantar Technology Inc. 16[. 
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Figure 1: Experimental setup used to measure the position resolu¬ 
tion of the Z-stack MCP-RA detector m- See text for details. 


Further details of the experimental setup can be found in 

i3- 

The entire assembly presented in Fig. [1] is housed in 
a vacuum chamber that is evacuated to a pressure of 4 
X 10“® torr. The microchannel plates were biased to a 
voltage of +3139 V using a ISEG NHQ224M low-noise, 
high voltage power supply (HVPS). The RA was biased 
to +3284 V, also using a ISEG NHQ224M HVPS. The 
secondary electron emission foil and photomultiplier tube 
were biased to voltages of -300 V and -1800 V using HK 
5900 and Bertan 362 HVPS respectively. Signals detected 
at each corner of the resistive anode were ainplified by a 
high quality charge sensitive amplifier (CSA) [l 8 | operated 
in vacuum. The four GSAs are situated approximately 13 
cm from the MGP-RA to mimimize cable capacitance. The 
output signals from the GSAs are coupled to a 250 MS/s 
digitizer (Caen DT5720B) 0 . Readout of the digitizer is 
triggered by the coincidence of a fast signal extracted from 
the back of the MCP detector and the PMT. The digitizer 
is readout by a standard PC and waveforms are recorded 
for subsequent analysis. Eor the analysis subsequently de¬ 
scribed the waveforms associated with a total of 260,000 
coincident triggers were recorded. 

Indicated in Eig. [5] is the reverse pincushion shape of 
the resistive anode (black) along with the circular outline 
of the 40 mm diameter MCP detector (blue). Superim¬ 
posed on the RA is the image of the ten slits provided by 
the stainless steel mask (green). The four corners of the 
MCP-RA are designated Qi, Q 2 , Q 3 , and Q 4 as evident in 
Pig. [5J Along with the relative position of the MCP-RA 
and mask, shown in Pig. [2]are the signals measured at two 
locations on the RA. One set of digitized traces (red) cor¬ 
respond to an electron cloud incident at the the bottom slit 
of the RA as indicated by the red square. The other set of 
traces (black) correspond to a position close to the center 
of the RA as indicated by the black square. The wave¬ 
forms associated with these two positions are markedly 
different. When the position signal arises from the center 
of the detector all the waveforms are essentially the same 
as is expected. However, when the signal originates at the 
bottom of the RA, while the two corners nearest the signal 
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Figure 2: (Color Online) Schematic of the RA (black) along with the 
outline of the 40 mm diameter MCP (blue) and slits in the stainless 
steel mask (green). Pulse shapes associated with the electron cloud 
incident on two locations of the RA are also shown. Slit 0 (sO) 
corresponds to the bottom slit and slit 9 (s9) to the topmost slit. 


exhibit a fast risetime followed by an exponential decay, 
the upper corners of the RA manifest a significantly slower 
risetime. This risetime, dictated by the RC of the resistive 
anode thus provides a measure of the particle’s position. 
Motivated by prior work which utilized the risetime of sig¬ 
nals in resistive silicon detectors to achieve position sen¬ 
sitivity in one dimension (^ . we elected to characterize 
the RA waveforms by their risetime. The risetime (RT) 
of each signal was defined as the time required for the sig¬ 
nal to go from 10% to 90% of its maximum value. As the 
slits in the mask are oriented to probe the Y dimension 
of the RA, for the remainder of this work we focus on the 
position information in that dimension. 


3. Signal Risetime analysis 


Signals obtained from the CAEN digitizer are processed 
through a series of mathematical operations using a stan¬ 
dard C++ code calling ROOT 21| libraries. In the present 
investigation the sampling resolution is 4 ns since the dig¬ 
itizer used has a sampling frequency of 250 MS/s. The 
obtained signals are corrected for the DC offset and gains 
for different channels for CSA’s on an event-by-event basis. 

Presented in Eig. |3] is the correlation between the rise¬ 
time observed for signals at the two bottom corners if the 
RA namely Q 3 and Q 4 . Given the exponential depen¬ 
dence of the signal amplitude on the RC of the resistive 
anode, the correlation is examined on a logarithmic scale. 
Individual slits are clearly evident in Eig. [3] with low num¬ 
bered slits exhibiting shorter risetimes and higher num¬ 
bered slits associated with longer risetimes. An overall 
linear dependence between log(RTQ 3 ) and log(RTQ 4 ) is 
observed as one moves from the bottom of the RA (sO) 
to the top (s9). Interestingly, a large jump in risetime is 
observed between slit 5 (s5) and slit 6 (s 6 ), indicating a 
high sensitivity of the risetime to position in the center 
of the detector. Moreover, these two slits, in contrast to 
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Figure 3: (Color online) Correlation between the risetime extracted 
from CSA output signals for the corners Q3 and Q4. The numbering 
of the slits is indicated in the figure. 


the other slits manifest a broad range of risetimes in at 
least one of the two risetimes. Electrons associated with 
these slits often do not show a strong correlation between 
the risetime of Q 3 and the risetime of Q 4 as indicated by 
the horizontal and vertical bands in the figure. Both the 
large sensitivity to position of the risetime in this region 
of the detector and the large dispersion in risetime are 
likely due to the fact that the charge cloud incident in this 
region experiences near equal resistance in all directions 
and a relatively small gradient in all directions. For all 
but these center two slits however the risetime measured 
for the two corners exhibits a clear anti-correlation. This 
anti-correlation arises from the spatial extent of each slit 
in the x dimension. Motivated by the overall positive cor¬ 
relation between log(RTQ 3 ) and log(RTQ 4 ) we construct 
the quantity log(RTQ 3 ) -|- log(RTQ 4 ) corresponding to a 
line of unity slope in Fig. [3] It is clearly evident that the 
spacing of the slits along this line is not uniform. 

In FigHlthe dependence of the summed risetime, log(RTQ 3 ) 
+ log(RTQ 4 ), on the position in the Y dimension is ex¬ 
plored. The position in the Y dimension is determined by 
using the charge division method proposed by [ 1 , @| . 


^position — (^3 “b /Qtotal ( 1 ) 

Qtotal = Ql + Q2 + Q3 + Qa ( 2 ) 

To obtain the charge measured at each corner of the 
MCP-RA we have utilized digital filtering techniques to 
provide signal conditioning. We initially used a gaussian 
filter on the digitized signal to integrate the signal from 
each CSA resulting in a gaussian-like pulse. The ampli¬ 


tude of this pulse is related to the charge measured at the 
corner. To efficiently realize the gaussian filter a recursive 
algorithm was employed 22, 2^ . Use of an integration and 
differentiation time constant of 500 ns was determined to 
result in the best position resolution. The resolution ob¬ 
tained is comparable to the resolution of 170 /im obtained 
for this experimental setup with analog electronics [13. 
To determine the sensitivity of our result to the filtering 
technique chosen we also used a trapezoidal filter. Fur¬ 
ther details on the signal processing can be found in the 
Appendix. 

Evident in Fig. |4] is a clear correlation between the 
summed risetime and the Yposiuon of the electron cloud. 
Since Yposition=0 is associated with the top of the RA, 
as the value of YposiUon increases one observes a general 
decrease in the summed risetime. This trend is under¬ 
standable since as the position of the electron cloud moves 
closer to the bottom of the RA the risetime decreases as 
initially evident in Fig. [T] The correlation observed for 
the upper half of the RA, YposiUon <0.45 is also observed 
for the lower half of the RA, Yposition> 0.45. The fact 
that the observed locus is not centered on Yposition=0-5 is 
likely due to the uncertainty in the positioning of the mask 
relative to the RA. Situated along the locus of points in 
Fig. m that is the main feature of the spectrum, are a 
series of peaks corresponding to the slits in the mask. An 
enlarged region centered on slit 3 is displayed in the inset 
of the figure. The relationship between summed risetime 
and Yposition exhibited by the peaks in 2 ] provides useful 
information. It describes how these two quantities should 
be related for a given position of the electron cloud on the 
RA. Selection of the data in these peaks corresponds to 
approximately 25 % of the data. Aside from the strong lo¬ 
cus with peaks visible in Fig. |4]one also observes data that 
does not lie on this locus. The majority of this data is as¬ 
sociated with horizontal lines which originate from one of 
the ten peaks. Examination of the pulse shapes associated 
with this data clearly establishes that much of this data 
is associated with the pileup of two signals. Such a pileup 
disturbs the baseline for the signal distorting the measured 
charge while minimally disturbing the signal risetime. 

To investigate if using the joint information from the 
summed risetime together with the charge division method 
provides an improvement in the position resolution we se¬ 
lected the ten peaks visible in Fig. 21 As the resulting 
distributions are reasonably well described by a gaussians, 
they were fit with this functional form. As the distance 
between the slits corresponding to each of the peaks is 
well known the position spectrum was calibrated and the 
position resolution extracted. The resulting resolution is 
presented in Fig. [SI 

To begin we examine the measured resolution where 
only the charge division technique is utilized for the digi¬ 
tized signals. The result for this case which does not make 
use of any risetime information is presented as the cross 
symbol (green). The resolution obtained from this analysis 
with the trapezoidal filter is typically 190 ^m (FWHM). 
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Figure 4: (Color online) Dependence of the summed risetime of 

Q3 and Q4 on the Y position of the electron cloud. The peaks 
corresponding to slit 0 (sO) and slit 9 (s9) are indicated. The peak 
associated with slit 3 is presented in the inset. 


Yposition (mm) 



-20 

-10 

0 

10 

20 


250 

■ □ 

_ + 

Analog 

Trap(CD) 



- 



- A 

Trap(CD+RT) 



- 



• 

Gaus(CD+RT) 




E 

200 



-i- i 

i ^ 






I 

T 

- 







_ 

H 

150 



9 










LL 






- 


100 


# ^ i 

i ♦ 

• ♦ 

* 

1 , , 

t - 



0 

2 

4 

6 

8 


slit number 


Figure 5: (Color online) Comparison of the measured resolution 
as a function of slit (position) using the correlation between risetime 
and position with the charge division approach alone. The resolution 
obtained using analog electronics is shown for reference. 


Use of the gaussian filter (not shown) gives a comparable 
result. As the use of these filters mimics the use of an 
analog shaping amplifier, also shown for comparison is the 
resolution obtained from the use of analog electronics. Al¬ 
though the resolution obtained with the analog electronics 
is Ri20 /rm lower than that obtained with the trapezoidal 
filter it is not a dramatic reduction. Moreover, using the 
analog electronics a reasonable resolution could only be 


obtained for the four central slits [17|. In contrast, us¬ 


ing the digital signal processing approach a relatively uni¬ 
form resolution is obtained across the entire detector. The 
somewhat lower value observed for the two edge slits is due 
to the slight deviation of the distributions from gaussians. 

The measured resolution obtained by the joint use of 
summed risetime along with charge division is depicted by 
the red circles (gaussian filter) and blue triangles (trape¬ 
zoidal filter) respectively. Both filters provide a measured 
resolution of ~90 p,m across entire detector. This result is 
a significant improvement over the use of the charge divi¬ 
sion approach alone. No systematic advantage is observed 
for one filter as compared to the other filter. Averaging 
the measured resolution of the central eight slits for either 
the trapezoidal or gaussian filters, one obtains an average 
measured resolution of 92 ^m. Accounting for the finite 
slit width of 100 ^m this measured resolution corresponds 
to an intrinsic resolution of 64 /rm . This result is com¬ 
petitive with the resolution obtained using more complex 
MCP arrangements and a retarding potential 11, MEl- 


4. Conclusions 

Resistive anode (RA) MCP detectors are widely used 
to provide position resolution in the detection of electrons, 
photons, or ions. A charge division approach from the four 
corners of the RA is typically utilized to extract the cen¬ 
troid of the electron cloud emanating from the MCP. The 
signals arriving at the four corners of the RA manifest a 
broad range of risetimes. For a 40 mm diameter MCP-RA 
detector, risetimes range from < l/xs to approximately 40 
/iS. By digitizing the signals from the MCP-RA with a 
high speed digitizer, digital signal processing techniques 
were utilized to extract both the charge collected at each 
corner as well as the risetime of the four signals for each 
incident electron. Implementation of the charge division 
approach resulted in a resolution of Ril90 ^m (FWHM), 
comparable to that obtained with analog electronics. Ex¬ 
amination of the digitized signals reveals a clear correlation 
between the risetime of the signals and the position of the 
electron cloud. By using this correlation between risetime 
and position together with the charge division method a 
measured position resolution of ~90 ^m (FWHM) was 
achieved, corresponding to an intrinsic resolution of 64 /rm 
(FWHM). This result represents a significant improvement 
in the position resolution obtained. It was also observed 
that in the central region of the detector the risetime is 
extremely sensitive to the position of the electron cloud. 
While this latter observation requires further investigation 
for full characterization, it presents intriguing possibilities 
for enhanced imaging capabilities. 

A key question in image processing is the speed with 
which signals can be processed. We therefore assessed the 
rate at which signals could be processed through our imple- 
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mentation of the digital filters described. Using a 2.8 GHz 
i686 PC we measured a rate of approximately 350 coinci¬ 
dent events/s for the trapezoidal filter and 200 events/sec 
for the gaussian filter. In the near term this processing 
rate could be increased by parallelization to make use of 
multiple cores on the computer. A longer term goal would 
involve implementation of FPGA processing at the detec¬ 
tor level to both achieve a high processing rate and reduce 
the volume of data acquired. 

5. Acknowledgments 

We gratefully acknowledge the technical support pro¬ 
vided by the personnel in the Mechanical Instrument Ser¬ 
vices and Electronic Instrument Services at the Depart¬ 
ment of Chemistry, Indiana University. This research is 
based upon work supported by the U.S. Department of 
Energy, National Nuclear Security Administration under 
Award Number de-na0002012 


relation [^. We utilized the pole-zero recursive relation 
given by the following equation 


HPLrrl^] 


= HPUi] + PZCxCSA^[t-l] (7) 


where PZC is a pole-zero correction factor. The value of 
PZC can be chosen so as to correct for the overshoot. The 
second term provides the amplified (attenuated) exponen¬ 
tial tail of the input pulse which corrects against the pole 
present in the exponential part of The pole-zero 

corrected output of the high-pass filter is then transformed 
by a four-pole, low-pass digital filter, which is realized by 
the following recursive relations 221 


GSHi\i] = + blpGSW[i - 1] (8) 

+blpGSH^[i - 2] -b blpGSW[i - 3] -b b'ipGSH^li - 4] 


6. Appendix: Signal Processing 


Standard digital signal processing techniques are used 
on the CSA signals to extract the position resolution using 
gaussian and trapezoidal filter algorithms 22, 2^ 
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To efficiently process the pulse shapes, recursive relations 
of both the hlters were employed. For the gaussian shaper, 
the input signal is subjected to a combination of high-pass 
and low-pass filters. Each sample, i of the channel 
CSA is subjected to a single-pole, hi gh-p ass digital filter 
using the following recursive relation 221 


where a°p, bpp, b\p, bpp, b\p are the coefficients of a 
low-pass filter kernel and GSH^ refers to the semi-gaussian 
shaper output obtained for the channel. The corre¬ 
sponding coefficients for the low-pass filter is given as 



(9) 

ilp = 

(10) 

blp = -6(e-i/^-)2 

(11) 

&ip = 4(e-i/^-)3 

(12) 

&PP = -(e-i/A.)4 

(13) 


HPLM = a%pCSA^[i] + a]jpGSA^[i - 1 ] 

+b\jpHPlut[i — 1] (3) 

where aPpp , Upp, and fepp are the coefficients of the high- 
pass hlter kernel, which can be calculated from the follow¬ 
ing equations 


“HP — (1 + e ^/"^^)/2 

(4) 

app = —(1-be ^/■^^)/2 

(5) 

b\,p = e-VAn 

(6) 


Ap refers to the number of samples being considered for 
the high-pass time constant. The factor decides 

the sample-to-sample decay for a given time constant. Equa¬ 
tion ([3]) mimics an analog differentiator circuit whose RC 
constant can be controlled by using the parameter Ap (in 
the units of number of samples). A time constant of 500 ns 
was used for all the channels of CSA. Due to finite length 
of the exponential tail of input signal, the high pass fil¬ 
ter output may overshoot the baseline of the pulse. This 
tendency can be corrected by using a pole-zero recursive 


where Ap is the low-pass time constant of the high-pass 
filter. Eor each corner of the resistive sheet, the measured 
charge is obtained by determining the height of the semi- 
gaussian shaper output. 

To ensure that the extracted resolution was not lim¬ 
ited by the particular filter chosen, we also implemented a 
trapezoidal digital filter to extract the charge of the input 
signals. The following recursive relations are used for the 

channel CSA 


= CSA^[i\ - CSA^[i - k\ 

(14) 


(15) 

p>\i\ = p^\i — 1] -b 

(16) 

r^[i] = p^[i] + 

(17) 

gJ [i] = g.? [i - 1] + [i] 

(18) 


Here k represents the number of samples in the rising edge 
while I is the total number of samples in the rising as well 
as in flat top region of the trapezoid shaped signal, m 2 is 
the gain of the filter, mi is the pole-zero correction factor, 
and is the output of the filter. The correction factor 


5 









depends on the decay time constant (r) of the preamplifier, 
given by the following equation [ 2 ^ 


™ _ ”^2 

“ eiTciklr) _ 1 


(19) 


where Tdk is the sample resolution of the digitizer. For 
shaping the CSA signals both the rising and flat top lengths 
are chosen to be 500 ns and the sample resolution is taken 
to be 4 ns. The gain of the shaper is taken to be unity, 
while the decay time constant is chosen to be 30 fis to 
match the decay constant of the CSA. Charge collected at 
each corner of the RA is calculated by taking the average 
height of the trapezoidal output. 
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